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ABSTRACT 

We discuss images of M87 from 3.6 to 160 /im obtained with Spitzer. As found previously, there 
is an excess in the far infrared over a simple power law interpolation from the radio to the resolved 
nonthermal features in the mid-infrared and optical. We show that this excess is most likely warm dust 
in the galaxy itself, and that the properties of this emission component are similar to the far infrared 
emission of normal giant elliptical galaxies. The new observations let us determine the spectrum of 
the jet and surrounding lobes of nonthermal emission. We find that even in the lobes the synchrotron 
break frequency is in the optical, probably requiring in situ particle acceleration not only in the jet 
but in the lobes as well. 

Subject headings: galaxies: active — galaxies: jets - galaxies: individual(M87) 



1. INTRODUCTION 

Fifty years ago, Baade (1956) discovered that the 
highly collimated jet in M87 (NGC 4486, 3C 274, Virgo 
A) glows in synchrotron emission. Since then, we have 
learned that such jets are often the most spectacular 
manifestation of an active galactic nucleus (AGN) and 
that the resulting beaming of energy strongly influences 
our view of the entire AGN phenomenon. Bec ause it 
is rel atively nearby (d=16Mpc, 78 pc arcsec -1 ; iTonrvl 
Il991|) . the jet in M87 remains one of the most readily 
studied of these structures. Much of our understanding 
of jets rests on the high- resolution, deep images of it now 
available from the radio through the X-ray. 

In the radio, an inner region (l'.5x2'.0) is embed- 
ded ina halo with a diameter up to ~15.'0 (e.g. 
[ Mills! Il952t iBaade fc Minkowskil Fl954t iRottmann et all 
1996). This region is composed of the nucleus, jet, 
North-pre cedin g (Np) l obe and South-following (Sf) lobe 
l)Turlandlll975llForsterl [l980'). High resolution radio im- 
ages of the inner region have revealed compl ex structure, 
with knots, filaments, e dges and rings (e.g. lOwen et alJ 
1198(1 iHines et all 11989th A radio outflow originating 
from the nucleus terminat es at the boundary of the ra- 
dio halo ijOwen et al]l200f)ft . In the optical, Hubble Space 
Telescope (HST) images reveal ap parent superlumina l 
motions for several knots in the jet ijBiretta et alJ ll999) . 
The images also indicate differences between the mor- 
phologies of the op tical and radio synchro tron emission 
ijSparks et ali ri996: Pe rlman et al]l2001aj) . M87 sits at 
the cen ter of the h ot X-ray atmosphere of the Virgo 
cluster l)Fabianll 1994(1 . with the gas temperature decreas- 
ing toward the central region. However, simple cooling 
flow models are not consistent with recent X-ray ob- 
servations, such as the low cooling rat e and the single 
temperature phase of intraclust er gas ijBohringer et alJ 
120011 iMolendi fc Pizzolatoll200il e.g.). In the subarcsec- 
ond Chandra X-ray image, two spectacular outflow-like 
structures extend from the nucleus east and southwest, 
respectively, possibly created by buoyant bubbles (e.g. 
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lYoung et al1l2002l iForrnan et alJEffll . 

Ground-based mid-IR observations have only detected 
the emission from the nucleus and the brightest knots 
of the je t, and are consi s tent w it h synchrotron emission 
mode ls (|Perlman et al.l l2001bt iWhvsong fc Antonucci 
2004). The Infrared Space Observatory (ISO) mid-IR 
imag es have marginally detected the lobe ijXilouris et alJ 
l2004j) . Other features, such as the individual filaments in 
the lobe, are only detected weak ly in the near-IR, bu t not 
yet in the mid-IR llStocke et aljll981[ iSmith et al.lll983t 
iNeumann et all 119951 iPerlman et1dlll2001a|) . ISO and 
the Infrared Astronomical Satellite (IRAS) observations 
of M87 reveal excess far-infrared emissi on above the syn- 
chrotron interpolat ion from the radio l)Haas et alJl2004t 
iXilouris et alJ l2()04). However, the available information 
in the infrared currently lags far behind the exquisite and 
very deep images available in the radio, optical, ultravio- 
let, and X-ray. To help remedy this situation, we present 
infrared images of M87 at 3.6, 4.5, 5.8, 8, 24, and 70/im 
and photometry at 16 Mm obtained w ith the Infrared 
Array Camera CIRAC: lFazio et alJpOOH) and the Mu lti- 
band Imaging Photometer fMIPS: lRieke et alJl200H) on 
Spitzer. 

2. DATA REDUCTION AND ANALYSIS 
2.1. HST image 

We retrieved archived HST data (PID 8048, PI J. 
Biretta) taken with the Wide Field and Planetary 
Camera 2 (WFPC2). The images were obtained on 
2001 November 17 through filters F450W, F606W and 
F814W. The data were processed through the Post Ob- 
servation Data Processing System ( PODPS) pipeline t o 
remove bias and flat-field artifacts l)Biretta et alj fl996). 
Individual exposures at each wavelength were combined 
to remove the cosmic-ray events. The outputs of the four 
chips of WFPC2 were combined and the final image was 
rotated to put north at the top and east to the left. 

To obtain the optical image of the jet and lobe, 
galaxy light was subtracted in IRAF as in lPerlman et alJ 
l)2001a|K We masked out the optical jet and the Sf lobe as 
indicated by the 5 GHz radio image. We also masked the 
region without data because of the smaller FOV of the 
PC chip compared to the WFC one. The task ELLIPSE 
in IRAF was used to fit isophotes of the galaxy image 
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and a model image was constructed with task BMODEL. 
The optical nonthermal emission was obtained by sub- 
tracting the model image from the original image with 
task IMCALC. We then masked out the globular clus- 
ters and chip joining region based on the residual image. 
The above processes were repeated several times until 
the final galaxy-subtracted image did not exhibit ring- 
like structure. In this final image, we replaced each pixel 
of the masked region except for the jet and lobe with 
a simulated value based on the average and scatter of 
residual counts in the non-masked region with the same 
distance to the galactic center. The final uncertainty in- 
cludes the fluctuation in the residual image, flat-fielding 
errors (~1%) and zero-point errors (~1%). 

As shown in the right column of Figure ^ the optical 
jet is well detected at all HST bands after subtracting 
the galaxy light. In the region of the radio Sf lobe, there 
is an elongated feature that becomes more prominent at 
the longer wavelengths. This feature has morphology, 
location and position angle similar t o the radio compo - 
nent 6 seen in the 6 cm image (See Hincs et al. 1989). 
The optical spectrum of this feature has a ~ 2 (where 
f„ oc v~ a ), steeper th an that of the optical jet (a ~1, See 
iPerlman et al.l2001a|) but flatter than that of galaxy light 
(~ 2.5). The non-thermal spectrum of this feature indi- 
cates it is the optical component of the Sf lobe. The op- 
tical counterpart of filament 9 has a lso been detected by 
previous ground-ba sed observations l)Stiavelli et alll992t 
ISmrks et alllTflfll . 

2.2. SPITZER DATA 

Our MIPS observations (PID 82, PI G. Rieke) were 
made with the standard small field photometry mode 
and were reduced with the MIPS i nstrument team Dat a 
Analysis Tool (DAT) version 3.03 {Gordon et al.ll2005ft . 
At 70 /im, extra processing steps beyond those in the 
DAT were applied to remove known transient effects as- 
sociated with the 70 /im detectors to achieve the best 
possible sensitivity. Detector-dependent structures were 
removed by subtracting column averages from each ex- 
posure with the source region masked. In addition, a 
pixel-dependent time filter was applied (with the source 
region masked) to remove small pixel-dependent resid- 
uals. For all the bands, the background was subtracted 
using the sky level at the edge of each field. The resultant 
images were resampled by a factor of four and rotated to 
put north at the top and east to the left. The final MIPS 
images have fields of view (FOVs) of 7f5x7f5, 6f5x3.'0 
and 2.'5x6f0 at 24 /im, 70 /jm and 160 fim, respectively. 

To isolate the output of the jet and lobes at 24 and 
70/im, we subtracted the nuclear emission. We used the 
empirical point spread function (PSF) based on the ob- 
served images of stars and provided by the Spitzer Sci- 
ence Center (SSC) 3 . We aligned the PSF with the cen- 
troid of the image and determined the normalization by 
matching the average surface brightness within a half- 
FWHM radius. Although the PSF image has a FOV 
smaller than the image of M87, the contribution is negli- 
gible compared to the Poisson noise at regions outside of 
the PSF coverage. The original and nuclear-subtracted 
images are shown in Figure |3 The uncertainty in the 
nuclear-subtracted image was estimated by repeating the 

3 http:/ /ssc. spitzer. caltcch.edu/mips/psf.html 



subtraction with the PSF placed over a 4x4 grid around 
the centroid of the image. This uncertainty and the Pois- 
son noise were added quadratically to create the final 
noise image. 

IRAC data were taken from the Spitzer archive (PID 
3228, PI W. Forman). We used the post-basic-calibrated- 
data (post-BCD) images at a pixel scale of l'/2 pixel -1 . 
The resultant images were also rotated to put north at 
the top and east to the left. 

To obtain the non-thermal emission of the jet and lobes 
at the IRAC bands, the stellar emission was subtracted 
similarly to the HST image. Figure 01 shows the original 
image and residual image of M87 at the four IRAC bands. 
For some of the faint low-surface-brightness features, we 
found it impossible to extract reliable measurements at 
3.6 and 4.5/im because of the residuals from subtraction 
of the stellar emission. With this exception, we obtained 
reliable measurements for two positions along the jet and 
four in the lobes at all the Spitzer bands out to 24/im, 
a measurement of the jet and Sf lobe at 70/im, and an 
integrated measurement of the nuclear region and sur- 
rounding galaxy at 160/im. 

Because of the complexity of the images, we could not 
apply standard aperture corrections to our photometry. 
Instead, we used a variety of approaches optimized for 
each situation. The nuclear fluxes were measured by 
PSF fitting and calibrated against standard stars mea- 
sured the same way. To determine spectral slopes, we 
convolved all the images to the resolution of the low- 
est resolution image and then used the same photome- 
try approach on all. To obtain flux densities for the jet 
knots, we used the _ffSTF814W image as a template. We 
first did photometry in the selected aperture on the full- 
resolution image, and then repeated the photometry on 
the image convolved with the 24/im PSF. A comparison 
of the results yielded the appropriate aperture correction. 
We quote all the IRAC photometry at this resolution also 
for consistency. However, the nonthermal lobes are not 
well enough detected to use the F814W image in this 
way, and they are more extended than the knots, so we 
report photometry for them with no aperture correction. 
For IRAC measurements of the full M87 galaxy, we ap- 
plied extended source corrections in accordance with the 
procedure recommended by the SSC. MIPS photometry 
of the whole galaxy was corrected in accordance with 
smoothed STinyTim models of the image, which have 
been shown to fit the observed PSFs well (Engelbracht 
et al. in prep., Gordon et al. in prep., Stansberry et 
al. in prep.). Color corrections are also applied to MIPS 
photometry for a power law spectrum with an index of 
1.0. 

The measurements are summarized in Tables 1-3. 
Notes to the tables give details regarding the aperture, 
extended source and color corrections. 

2.3. Radio Data 

A VLA 5GHz radio map of M87 w as obtained at a res - 
olution of -0.4" on 1986 April 25 bv Mines et all (|1989|) . 
who describe the observation strategy and data reduc- 
tion. 

2.4. Image Convolution 

The images must be convolved to the same reso- 
lution as MIPS at 24/im to construct the full spec- 
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tral energy distributions (SEDs). Therefore, we con- 
volved the HST and the 6 cm radio data with the 
MIPS PSF. For the 6 cm image, we subtracted the 
nuclear emission by masking the central l'.'6xl'/6 re- 
gion before convolution. The IRAC images were con- 
volved with kernels that transform the observed IRAC 
PSFs to the MIPS 24 fim PSF (Gordon et al., in prep.) 
The noise images were calculated by error propaga- 
tion. The centroids of all images were aligned to the 
center of the 6 cm data at RA(J2000)=12h30m49.42s, 
DEC(J2000)=12d23m27.97s. 

3. RESULTS 

3.1. Infrared Image 

Figure H shows the radio, MIPS 24 /im, IRAC 8 /im 
and optical image of M87 with nuclear subtraction at the 
resolution of MIPS at 24 /im. The stellar emission was 
also subtracted for the IRAC and optical images. It is 
obvious that the jet, Np and Sf lobe are all detected at 

8 /im and 24 /im, with morphology similar to the radio 
emission. The infrared jet extends from the nucleus along 
a straight line until it bends toward the southeast near 
the boundary of the Np lobe. The infrared jet peaks 
at the position of knot-A (although not all knots are 
resolved), which is the brightest knot in the radio and 
optical. The brightness profile across the jet indicates the 
width of the jet is not resolved at 24 /im but marginally 
resolved at 8 /im, consistent with the width of the jet of 
around 1" at knot A in the optical image. As illustrated 
by contours in Figure 0] unlike in the radio and optical 
image where only one flux peak is present in the Sf lobe, 
there are two flux peaks near the bounda ry of the Sf 
lobe i n the infrared. The radio filaments (See lHines et alJ 
1989) corresponding to these two peaks are the filaments 

9 and rj. 

Besides the infrared jet and lobes, an infrared halo is 
clearly visible at 24 /im. The subplot in Figure [21 shows 
the surface brightness of the emission at 24 /im after 
masking the central lobe r egion defined by th e lowest 
level contour in Figure 4 of IHines et al.l l|1989f) . Unlike 
the lobe, it has a spherical structure as shown in Fig- 
ure [21 with flux decreasing with increasing distance from 
the nucleus. Moreover, although the infrared halo is on 
the same scale as the radio inner region (l'.5x2'.0), the 
morphologies are quite different since the radio inner re- 
gion is elongated and aligned with the radio lobe (See 
lOwen et al.ll20fl(T) . The discrepancy in the morphology 
implies that the emission of the infrared halo arises from 
a different mechanism. Figure [S] shows the map of the 
spectral index Q^q^. where the photometry at each grid 
point at 5 GHz and 24 /im was carried out on the im- 
age at the resolution of MIPS at 24 /tm and with the 
nucleus subtracted. The halo has ou,qh z flatter than the 
nonthermal emission region, again suggesting a different 
emission mechanism. 

As shown in Figure|21 at 70 /im, the jet and Np lobe are 
detected but are blended; the Sf lobe is also detected. At 
160 /im, the jet and lobes are all blended at the Spitzer 
resolution. The infrared halo is not detected at either 70 
/im or 160 /tm. 

3.2. Infrared Emission from the Host Galaxy 



Table ^ shows photometry of M87 within a radius of 
30", including emission from the nucleus, jet, and most 
of the Np and Sf lobes. The photometry for the Two 
Micron All Sky Survey (2MASS) was carried out on the 
original images without nuclear subtraction. The ISO- 
CAM photometry is from lXilouris et all |200l) ; the ISO- 
CAM measurements at shorter wavelengths have been 
supplanted by the IRAC data, which has a substantially 
improved ratio of signal to noise. We also include IRAS 
Faint Source Catalog (FSC) data. Since it was extracted 
using a point source filter, the effective apertures should 
be close to the 60" used for the other data sources. 

To interpret these data, we fitted a st ellar SEP in the 
near infrared, taking galaxy colors from I Johnson! l)1966f) 
(because this reference consisten tly integrat e s JKL - we 
applied color corrections as in ICarpenterl <|2001() and 
Rieke et al., in prep.). The SED was continued to the 
longer IR AC bands using the m easurements of individual 
K stars in Re ach et alJ (|2005j) . In normalizing the stel- 
lar spectrum, we emphasized the IRAC 3.6 and 4.5/im 
points because systematic errors within the IRAC mea- 
surements should be smaller than between them and data 
from other sources. Specifically, the photometry in all 
four IRAC bands was conducted on the SSC post-BCD 
images, using the same apertures and sky reference areas. 

We fitted two power laws to the nonthermal spectrum, 
because the n ucleus is self-absorbed at freque ncies below 
~2 GHz (e.g. iCharlesworth fc Spenceilll982j) and hence 
has a different spectrum from the other components. The 
n uclear component is fitt ed to the 10.8/im measurement 
of iPerlman et alJ l)2001b|h our measurements at 24 and 
70/tm, and those reported bv lHaas et alJ l)2004|) at 450 
and 850/xm. It has a — 1.15 and a flux density of 45mJy 
at 24/im. The remaining nonthermal emission remains 
optically thin into the radio, down to 1.4GHz. It has been 
fitted with a power law with a — 0.92 and a flux density 
of 50mJy at 24/im. Because the break frequency for the 
majority of the nonthermal emission is in the optical or 
UV (Section 3.3.3), we continue the power laws toward 
shorter wavelengths. 

Table^shows the aperture photometry, the stellar and 
nonthermal spectral components, and the residual flux 
densities when the modeled SEDs are subtracted from 
the measurements. This information is also illustrated in 
FigureEl There is a well detected excess above the model 
from 6 through 70/im, which is particularly strongly in- 
dicated at 24 and 70/im. However, the measured val- 
ues for the nonthermal emission (nucleus, jet, and lobes) 
at these two wavelengths fall very close to the model, 
which is dominated by the nonthermal emission. There- 
fore, it appears that the excess is not associated with 
the resolved nonthermal sources. The possibility that 
this emission is associated with the extended nonther- 
mal radio emission can be tested from Figure 5, which 
shows that the spectrum between 5GHz and 24/im be- 
comes flatter with increasing distance from the nucleus 
and jet. This behavior is not expected for extended non- 
thermal emission, which should either retain the same 
slope or have a steeper one due to energy losses in the 
s ynchrotron electron spectrum. 

Spa rks et al.l (|1903f) discuss filaments visible in Ha. It 
is plausible that dust in these filaments contributes some 
of the infrared emission not associated directly with the 
nonthermal source components. However, the morphol- 
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ogy of the 24/xm image (Figure 0) shows it to be more 
symmetric and more extended than the Ha filaments, 
which are concentrated in the region occupied by the 
jet and lobes. Therefore, most of the infrared emission 
must arise from a different component of the galaxy. The 
agreement with the pure e lliptical isophot es seen in the 
optical surface brightness ijLiu et alJl2005(l indicates the 
dust should be smoothly distributed in the host galaxy. 

The overall level of excess emission at 7 urn is at a 
normal level for giant elliptical galaxies (e.g. lTemill2004l 
lLeeuwi2004|) , where the dust is heated by the stellar light 
or collisions with electrons in the hot gas. The over- 
all SED of M87 from the near infrared to 24/xm is very 
similar to those of the low-activity (i.e., non-infrared- 
lu minous) brightest c luster galaxies (BCGs) discussed 
bv lEgami et alJ (|2006f> . The excesses detected weakly in 
the IRAC 5.8 and 8^m bands and in the IRAS FSC at 
12/im suggest the possibility of weak, extended aromatic 
emission, which is s een in a number of X -ray-emitting 
elliptical galaxies bv I Kaneda et all l|2005|) (but see also 
iBreeman et alJ ( f 2006'l ) . However, aromatic emission i s 
not seen in the nuclear spectrum ijBressan et al.ll200r^ . 
which instead shows silicate emission either associated 
with the AGN or with evolved stars. Therefore, if these 
features are present they must arise from dust distributed 
within the galaxy and perhaps heated as in other X-ray- 
emitting giant ellipticals. In any case, the characteristics 
of the mid- and far-infrared excess above the compact 
nonthermal emission from M87 appear to be consistent 
with expectations for a normal BCG and need have little 
to do with the nonthermal activity in the nucleus and 
jet. 

3.3. Behavior of the Nonthermal Sources 
3.3.1. SEDs of the Nonthermal Sources 

We now consider the extranuclear nonthermal emis- 
sion. As shown in Figure we constructed the SED 
from the radio to optical wavelengths for two regions in 
the jet (knot-A/B and knot-C/G 4 , one region (labeled 
as Np-1) in the Np lobe, and three regions (Sf— 1, Sf— 9 
and Sf— rf) in the Sf lobe. Each region is circular with a 
radius of 5". The position of each region and the corre- 
sponding 24 fim flux density are listed in Table 

Figure shows the SEDs of the different regions. The 
stellar emission is modeled with a blackbody spectrum 
fitted to the total IRAC 3.6 fim and total optical emis- 
sion within each region. Table lists the radio- infrared 
spectral index ot 2 ^" 1 between 5 GHz and 24//m after sub- 
tracting the stellar contribution and the radio-optical in- 
dex a^' between 5GHz and the three HST bands. In the 
jet and lobe, it seems that the regions with lower infrared 
surface brightness nave a steeper a R . a R is also 
flatter than a^', as expected for synchrotron radiation 
with a higher rate of electron energy loss at the higher en- 
ergies. The larger steepening Aa = a^' — a 2 R Aim in the 
lobe (-0.03-0.1) than the jet (-0.01-0.03) reveals that 
the high-energy cutoff of electron energy in the lobe oc- 
curs at a lower frequency, as demonstrated by the SED 
fits below. 

4 we use composite designations to emphasize that the knots are 
blended at our resolution. 



3.3.2. Minimum Pressure Analysis 

Assuming the equipartition condition that the parti- 
cle energy is equal to the magnetic energy, the mini- 
mum magnetic field, energy, and pressure of the emit- 
ting plasma can be derived from the observ ed SED, the 
emitt ing volume and the filling factor cf> llPacholczvld 
1970). The standard minimum pressure analysis assumes 
a power law distribution of electron energy; the depen- 
dence of the physical parameters on the high frequency 
cut-off is negligible if the spectral index a > 0.5, as for 
M87. This implies that the synchrotron loss at high fre- 
quencies does not change the results for M87. 

During the calculation, the ratio of proton energy to 
electron energy is assumed to be unity. The SED is as- 
sumed to be a power law with a spectral index equal 
to a^ 4/im between 10 s Hz and the break frequency A 
low frequency cut-off of 10 8 Hz was adopted, as the SED 
of jet and lobe starts to deviate from a powe r law at 1 
GHz (|Feltenlll96cl iMeisenheimer et all 11996). One or- 
der of magnitude variation in the low frequency cutoff 
results in a factor of —1.2 change in the magnetic field 
and a factor of —1.4 change in the energy and pressure. 
The volume is taken to be spherical with a radius of 5 
arcsec and the filling factor is assumed to be 0.1. A fac- 
tor of 10 variation in the filling factor corresponds to a 
factor of —2 in the magnetic field, energy, and pressure. 
The filling factor is hard to determine as it depends on 
the true three dimensional structure of the synchrotron 
plasma. The value of 0.1 was adopted to balance among 
various alternatives, such as a whole filled volume or a 
thin bo undary layer of th e volume as proposed by some 
studies l)Owen et al.lll989l) . 

The last three columns of Table list the minimum 
magnetic field, energy, and p r essure computed with for- 
mulae given by iBurns et alJ lll979ft Our results agree 
well with those of lEilek et alJ l|2003). Independent sup- 
port for the derived magnetic field strengths comes from 
the inverse Compton fluxes measured in Tev gamma rays. 
IStawarz et all pOOl analyze the HESS and HEGRA 
data to derive a lower limit of 300^G. However, the Tev 
gamma ray limit cannot be applied individually to the 
various regions we have isolated in M87. 

3.3.3. Model of the Synchrotron Emission 

We employ a quantative model to extract the un- 
derlying physical parameters regulati ng the SED. The 
Kardashev-Pacho lczyk (KP) model ijKardashevI Il962t 
lPacholcz vklll97 0|) has been shown to fit the jet behavior 
well l)Perlman et al.l2001a[) . It describes the synchrotron 
emission of an ensemble of electrons ejected t seconds 
ago. Assuming the initial energy distribution of the elec- 
trons is a power law N(Eo,oj, 0) = NqEq 7 where ui is 
the pitch angle, the new energy distribution of electrons 
at time t due to synchrotron loss is: 

N(E, uj, t) = NoE^il - c 2 B 2 sin 2 {uj)Et)~<- 2 (1) 

where c? = 2.37 x 10~ 3 (in cgs units) as defined in 
lPacholczvldl(T970]l . Integrating over the pitch angle, the 
intensity of the synchrotron emission from these electrons 
at frequency v is: 

I v = 2ir S c 3 c { r 1)/2 N B^+ 1 V 2 4 1 --< )/2 B(i T , 7 ) (2) 

where s is the extension of the source along the line of 
sight, Ci=6.27xl0 18 (in cgs units), c 3 =1.87x 10~ 23 (in 
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cgs units), B is the magnetic field, z>t = c\/{c^B 3 t 2 ). 
and xt = v/vt- Therefore, for a given magnetic field, the 
shape of the synchrotron radiation spectrum is described 
by three physical parameters: the break frequency ( i>t 
), the initial power-law index of the electron energy dis- 
tribution (7) and the normalization (flux at a given fre- 
quency) . 

The dotted line in Figure[7|shows the SED model fitted 
to the photometry excluding upper limits. The results 
from the SED fits for regions in the jet and lobe are given 
in Table 01 As a check of our calculations, the break fre- 
quency in knot- A/B agrees well with previous work; for 
example, it is intermediate between the values for the 
indivi dual knots A and B obtained by IWaters fc Zepfl 
( 2005). Our values for the knots are > 1 - 5 x 10 15 Hz. 
In the Np and Sf lobes, the turnover of the spectrum oc- 
curs in the optical with a break frequency at least a fac- 
tor of five lower than for knots A/B, except for the Sf— 9 
regi on, which has a com parable break frequency to the 
jet. iStiavelli et alJ l)1997fl obtained a cutoff frequency of 
4.3xl0 14 Hz for the feature 9, almost one order of magni- 
tude lower than our value. The discrepancy arises mainly 
because of the difference in the optical photometry of the 
feature 9. We believe that the high resolution and high 
S/N HST images should minimize the uncertainties asso- 
ciated with the stellar subtraction. Nonetheless, accurate 
determination of the break frequencies requires photome- 
try in the UV, which is not used in either study. However, 
our conclusions in this paper depend only on the break 
frequencies lying at optical or shorter wavelengths, not 
on the exact frequencies involved. 

The lowerlimit of the current maximum electron energy 
ljUchivama et a l. 2006) derived from the break frequency 

E max = 0.28(^) a5 (i?/10- 4 )- a5 (Vl0 14 ) - 5 Tev 

(3) 

is another physical parameter constraining the accelera- 
tion model, where z is the redshift and 8 is the Doppler 
factor. For z=0 and 8=1, as shown in Tabled] a max- 
imum electron energy with Lorentz factor greater than 
10 6 prevails in both the jet and lobe regions. 

3.3.4. In Situ Reacceleration in the Lobe 

The Spitzer data provide improved constrains on the 
synchrotron spectrum in the Np and Sf lobes, show- 
ing that the break frequency occurs in the optical or in 
the UV for Sf— 9. The c orresponding sync hrotron age 
(t = (ci/(c%B 3 i/ t )) - 5 sees: lPacholczvklll97'o|) for the two 
jet regions and Sf— 9 is around 500 yrs while electrons 
in the Np and Sf lobes are ~lxl0 3 yrs old. The largest 
distance that electrons can travel in 1000 yrs must be 
smaller than ct = 300pc (4") where c is the speed of light. 
Therefore, in situ reacceleration of electr ons must occur 
in the jet (as found previously by, e.g.. iPerlman et al.l 
I2001al lEilek et all 120031: IWaters fc Zepl 120051) and in 
both the Np and Sf lobes. iHines et all l|1989|) showed 
how this process could explain some additional features 
in the lobes, such as the bright radio filaments. 

As shown in Tabled the individual regions in the lobes 
have a minimum electron energy around 10 54 erg and 
thus a conservative estimate of the electron energy in an 
entire lobe is around 10 55 erg, considering its probable 
total volume. Given a synchrotron age of 10 3 yrs, the 



minimum input energy rate into an entire lobe is around 
3x 10 44 erg s _1 . Assume for the moment that this energy 
is supplied by supernovae. If each supernova releases an 
energy of 10 52 erg, the explosion rate accounting for the 
observed input energy rate is 1 yr , much higher than 
observed in normal galaxies. The gravitational poten- 
tial as implied by the X-ray luminosity of thermal gas 
is much lower compared to this reacceleration rate. The 
most plausible energy source is the jet, especially for the 
Np lobe where the jet is v isible. The minimum jet power 
is around 3xl0 44 erg s -1 ijOwen et al.ll2000|) . If the elec- 
trons in the lobe are reaccelerated in situ by a mechanism 
driven by the jet, most of the jet power must be converted 
into the kinematics of electrons in the lobe. 

Radial outflow connecting the lobe and outer halo 
is frequently interpr eted as a buoyan t bubble origi- 
nally inflated by jet ijOwen et alJl2000|) : the radio lobe 
may be i nflating the bubble befo r e detaching from the 
nucleu^jGjjll^^J^rttigy^B 119731 : iBicknell fc Begelmanl 
1996: iChurazov et alJ l20011. The rising bubble captures 
the hot X-ray gas from the ambient m edium, creating the 
X-ray arc as seen in the X-ray image l|Young et al.l l2002: 
iForman et al.ll2006h . In general, the total energy of jet 
power is converted to the internal energy of the lobe o r 
to the work to inflate the lobe (e.g. IBicknell et aLlll997|) . 
The relative fraction may be determined by the expan- 
sion process. The result of this study, that most of the 
jet power is locked into the kinetic energy of the elec- 
trons in the lobe, implies that only a small fraction of 
the jet power is used to inflate a cavity (radio lobe) in 
the thermal gas atmosphere. This conclusion appears to 
be consistent with the X-ray observation that the energy 
for i nflation is much lo wer than the minimum jet power 
fSee lYoung eXaDT 2002;i. The kinetic energy of the high- 
energy electrons should be carried away by radiation due 
to the short synchrotron lifetime, while the radio plasma 
could live long enough to inflate the lobe if synchrotron 
radiation is the only energy loss mechanism. As shown in 
Tableland Table 01 the total synchrotron luminosity of 
the lobe is much smaller than the energy rate for electron 
reacceleration. It is possible that most of the low-energy- 
electron kinetic energy is carried away by bubbles to the 
outer halo to heat the cooling flow. 

4. CONCLUSIONS 

We present Spitzer observations of the giant elliptical 
galaxy M87 in the Virgo cluster. The basic conclusions 
are: 

(1) The far infrared excess emission above the non- 
thermal power law is likely to be warm dust distributed 
throughout the galaxy. Its luminosity is similar to that 
observed in the far infrared for other, normal giant ellip- 
tical galaxies. 

(2) The infrared emission in the jet and two lobes is 
synchrotron emission with break frequencies in the UV 
(jet) or optical (lobes). The high break frequency in the 
lobes indicates that in situ reacceleration of electrons oc- 
curs there as well as in the jet. The minimum input 
energy rate to reaccelerate the electrons in the lobe is 
comparable to the jet power, implying that most of the 
jet power is converted to electron kinematics in the lobe, 
not used to inflate the lobe. 
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TABLE 1 

Photometry within the central 60" diameter region 



Wavelength 


Total Flux Density 


Stellar Component 


Nonthermal Component 


Residual 


(H 


(mJy) 


(mJy) 




(mjy) 


(mJy) 


1 oc; 


1 OA 1 _L1 OA 1 




5 






2.16 


1202±120 1 


1110 


9 




- 


o.O 


DUO 


DIG 


15 






4.5 


375 2 


365 


19 






5.8 


302 2 


225 


21 




53±15 6 


8 


189 2 


121 


33 




35±10 6 


12 


231±37 3 


55 


49 




127±40 


15 


106±21 4 


35 


62 




9±25 


24 


171.9±13.4 5 


11 


99 




59±15 


60 


394±51 3 


2 


251 




141±60 


70 


455.2±9.2 5 


1 


293 




161±15 


160 


581.5±10.0 5 





682 







Note. — 1 2MASS; 2 IRAC, fro m PIP 3228, PI W. F orman. We used extended source corrections of 0.934 at 3.6/im, 0.97 at 4.5/im, 0.86 at 5.8/Jin 
and 0.80 at 8/im; 3 IRAS FSC; 4 IXilouris et all 120041 : 5 MIPS, from PID 82, PI G. Riekc, aperture corrections of 1.09, 1.31 and 1.93 respectively 
at 24, 70 and 160 color corrections of 1.04, 1.1 and 1.04 respectively at 24, 70 and 160 (im; 6 Er rors are taken to be 5% of the total flux density 
in the band. 



TABLE 2 

Photometry at 70 fiM of the nucleus, jet and lobes 



Region 


RA 


DEC 


/70nm(Jy) 


Nucleus 


12 30 49.42 


12 23 27.97 


228 1 


Jet & Np lobe 


12 30 48.21 


12 23 33.9 


70 2 


Sf lobe 


12 30 50.49 


12 23 14.2 


57 2 



Note. — 1 From PSF fitting; 2 Aperture correction factor of 1.3 for a 30"-diameter aperture; Color correction factor of 1.1 for all three fluxes. 



TABLE 3 

Spectral properties of the nucleus, jet and lobes at the resolution of MIPS at 24/zm 



Region 




RA 




DEC 


/24^m 


fS/J,m 




/4.5/im 




24um 
a R 


opt 
a R 
















(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(erg) 




Nucleus 










50.7 1 


















knot-A/B 


12 


30 48.64 


12 


23 32.2 


32.1 2 


13.19 


10.29 


8.92 


8.36 


0.70 


0.71 


1.81E+ 


42 


Knot-C/G 


12 


30 48.12 


12 


23 35.4 


17.6 3 


5.97 


4.77 


3.84 


3.27 


0.74 


0.77 


6.88E+ 


41 


Np-1 


12 


30 47.88 


12 


23 24.5 


3.22 4 


0.81 


0.51 


0.54 


0.44 


0.85 


0.91 


1.03E+ 


41 


Sf-1 


12 


30 50.19 


12 


23 15.3 


3.45 4 


0.47 


0.57 






0.89 


> 0.99 


9.29E+ 


-10 


Sf-0 


12 


30 50.85 


12 


23 17.5 


2.93 4 


0.70 


0.63 






0.91 


0.94 


1.01E+ 


41 


Sf-T? 


12 


30 50.35 


12 


23 5.2 


2.70 4 


0.57 


0.46 






0.88 


> 0.98 


7.85E+ 


-10 



NOTE. — 1 Nuclear flux density obtained by fitting PSF to the surface brightness profile; 2 An aperture correction factor of 1.76; 3 An aperture 
correction factor of 1.40; 4 No aperture correction; Color correction factor of 1.04 at 24/^m. 
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TABLE 4 

The physical parameters from the spectrum fit 



Region 


aiNPUT 


Kt(Hz) 


tsyn (y r ) 


Ee max (0.5Mev) 


B m inp(G) 


EminP (erg) 




(dyn cm 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


knot-A/B 


0.70 


> 5E15 


< 5E2 


> 4E06 


1.3E-04 


3E54 


8E-10 




Knot-C/G 


0.73 


> 1.4E15 


< 8E2 


> 2E06 


1.2E-04 


2.5E54 


6E-10 




Np-1 


0.85 


~ 1.6E15 


~ 9E2 


~ 2E06 


1.0E-04 


1.6E54 


4E-10 




Sf-1 


0.90 


~ 6E14 


~ 1.4E3 


~ 1.3E06 


1.1E-04 


2E54 


5E-10 




Sf-0 


0.90 


> 2E15 


< 7E2 


> 2E06 


1.1E-04 


2E54 


5E-10 




Sf-7J 


0.88 


~ 6E14 


~ 1.5E3 


~ 1.4E06 


1.0E-04 


1.6E54 


4E-10 






Fig. 1 — 



Original HST images of M87 in the left column and galaxy-light subtracted image in the right column. 
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MIPS24/xm 159"X159" 




Radius (arcsec) 




I i I o.o \ ,, A 

5 10 15 20 25 30 
Radius (arcsec) 

Fig. 2. — Original MIPS images of M87 in the left column. MIPS images of M87 after subtracting the nucleus in the middle column. 
The radial profile of surface brightness is shown in the right column where the solid line is the profile of the original image and the dotted 
line is the profile of the PSF at the corresponding wavelength. The subplot at 24/tm shows the su rface brightness of the emission after 
masking the central lobe region defined by the lowest level contour in Figure 4 of Hincs ct al. 1 1989). 
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IRAC- 8^m (84"X84 M ) 




knot-A/^^' 



knot-C/G 




l 



IRAC- 5.8/xm (84"X84") 




IRAC- 4.5/xm (84"X84' 





IRAC- 3.6/xm (84"X84 I 




Fig. 3. — Original IRAC images in the left column and galaxy-light subtracted images in the right column. 
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VLA-6cm 68"X59" MIPS-24^m 68"X59 




IRAC-8/mn 68"X59" HST-F814w 68"X59" 




Fig. 4. — The radio, MIPS 24^m, IRAC 8/im and HST optical image of M87 after subtracting the nucleus at the resolution of the MIPS 
24fim image. The galaxy light is also subtracted for the IRAC and HST images. The contours (thin solid lines) arc superposed to illustrate 
two local maxima in the Sf lobe at the infrared wavelength while only one maximum is present in the radio and optical. 
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+ 12 23 54.5 



-0.59 0.54 0.56 0.55 0.58 0.61 0.67 0.73 0.71 0.62 0.56 



0.62 



+ 12 23 30.0 



o 

O 



+ 12 23 5.5 



0.67 



.0.82 



0.79 



"0.77 



0.67 



0.58 




-0.55 0.66 0.78 0.79 0.73 0.63 0.64 0.63 0.61 0.61 0.62 



+ 12 22 41.01 i — 

1 2 30 52.1 



12 30 50.0 12 30 47.9 

RA 



12 30 45.8 



Fig. 5. — The map of the spectral index c?^qh z - The measurements of flux density at 5GHz and 24 (im were made on the image with 
nuclear subtraction at the resolution of MIPS at 24 (im. The solid contour is the intensity distribution of the original 6 cm image. Contour 
intervals are [0.11, 0.43, 0.53, 0.64, 1.06, 5.32, 10.64, 173.86] mjy. 
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10 12 10 13 10 14 

Frequency (Hz) 

Fig. 6. — Comparison of large aperture photometry with models of the stellar and nonthermal SEDs. The light lines show the models 
of the power law nonthermal emission and the stellar output, and the heavy line is the sum of the two components. The open circles are 
the 60 arcsec diamet er aperture photometr y. The filled boxes show the flux densities for the nucleus, jet, and lobes from this work at 70 
and 24/^m, and from lPerlman et al. (2001b) (with a correction for low surface brightness components) at 10.8/^m. The diamonds show the 
residual fluxes after subtracting the model from the totals measured. At 70^tm, we label the points for total, nonthermal, and residual flux 
density. 
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Fig. 7. — SEDs of individual regions in the MIPS 24 (im image (see Figure[3J. The photometry of the nonthermal features is shown as 
diamonds, with downward pointing arrows for the upper limits. The lines are the fits of the KP synchrotron emission models. 



